Speckle interferometry 



> 



- - - Swapan K Saha 

^Zf^ . Indian Institute of Astrophysics 

Q^ ; Bangalore 560 034 

India 
e-mail: sks@iiap.ernet.in 

o ■ 
O 

oo _ 

CN , Abstract : We have presented the basic mathematical treatment of interferom- 

etry in the optical domain. Its applications in astronomical observations using 
both the single aperture, as well as the diluted apertures are described in de- 
\j^ • tail. We have also described about the shortcomings of this technique in the 

^^ i presence of Earth's atmosphere. A short descriptions of the atmospheric tur- 

^O ' bulence and its effect on the flat wavefront from a stellar source is given. The 

^— ^ , formation of speckle which acts as carrier of information is defined. Laboratory 

experiments with phase modulation screens, as well as the resultant intensity 
^>J , distributions due to point source are demonstrated. The experimental method 

to freeze the speckles, as well as data processing techniques for both Fourier 
modulus and Fourier phase are described. We have also discussed the technique 
of the aperture synthesis using non-redundant aperture masks at the pupil plane 
O I of the telescope, emphasizing set on the comparison with speckle interferome- 

-4—' ' try. The various methods of image restoration and their comparisons are also 

j^ I discussed. Finally, we have touched upon certain astrophysical problems which 

can be tackled with the newly developed speckle interferometer using the 2.34 
meter Vainu Bappu Telescope (VET) , situated at the Vainu Eappu Observatory 
(VEO), Kavalur, India. 
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1. Prologue 

When a person looks at the medium size optical telescope of 2 meter class, he 
may wonder about its large structure rather than the complexity involved in 
making a good mirror. Much more complexity is in store for an astronomer 
when he quests about the difficulties in obtaining diffraction limited images. 
Where is the hindrance? Is it that the aberrations persist in the telescope or 
the atmospheric turbulence creating the problem or both? 

Will they achieve the goal if a large telescope (diameter > 2 meter) is installed? 
Large telescope helps in gathering more optical energy, as well as in obtain- 
ing better angular resolution. The resolution increases with the diameter of 
the aperture. Owing to the diffraction phenomenon, the image of the point 
source (unresolved stars) cannot be smaller than a limit at the focal plane of 
the telescope. This phenomenon can be observed in ocean, when regular waves 
pass through an aperture. It is present in the sound waves, as well as in the 
electro-magnetic spectrum too starting from gamma rays to radio waves. 

Though interferometry at optical wavelengths in astronomy began more than a 
century and a quarter ago [1], but real progress has been made at radio wave- 
lengths in post war era. Development of long baseline and very long baseline 
interferometry (VLBI) , as well as usage of sophisticated image processing tech- 
niques have brought high dynamic range images with milliarcseconds (marcsec) 
resolution. 

The understanding of the effect of atmospheric turbulence on the structure of 



stellar images and of ways to overcome this degradation has opened a channel 
to diffraction limited observations with large telescope in the optical domain. 
Speckle interferomctric technique [2] is being used to decode the high angular 
resolution information. Diffraction limited information can also be obtained 
with a large telescope using other techniques viz., (i) pupil plane interferometry 
[3], (ii) differential speckle interferometry [4], (iii) phase closure technique [5], 
(iv) aperture synthesis using partial redundant, as well as non- redundant mask- 
ing method [6 - f5], (v) speckle spectroscopy [16], (vi) speckle polarimetry [17] 
etc. 

In the optical band, a large mirror oi (p >10 meter class with a high precision 
accuracy in figuring may not be possible to develop, therefore, the resolution 
is restricted with the size of the telescope. Introduction of long base line in- 
terferometry using diluted apertures became necessary and the improvements 
of instrument technology made it possible to obtain a very high resolution of 
the order of a few marcsec [18, 19]. Recent success in implementing phase- 
closure technique on three ground based telescopes in the optical domain [20] 
has produced new results with very high angular resolution. 

Developments of high resolution imaging have been going on in our institute 
over a decade. These were in the form of concentrating on theories of image 
analysis, as well as on conducting several experiments at the various telescopes 
and at the laboratory [21-26] and Recently, we have developed a basic speckle 
interferometer [27, 28] for use at the 2.34 meter VBT, VBO, Kavalur. 

Due to the paucity of funds, limitation of infrastructure facilities and the lack 
of support that includes man power etc., this project could not progress at a 
faster pace. Quite a bit of time had been spent on R & D during the execution 
of the project [21, 22, 24, 27, 29-31]. A few vital optical components had to 
be developed by us [27]. Several algorithms were developed simultaneously 
with bare facilities to carry out image analysis and processing [23, 26, 32, 33]. 
Nevertheless, I could complete the mammoth task successfully. [I remember a 
few words of a scientist, 'one has to have an infinite patience to undertake any 
new challenging project and bring it to a success'. I do realize now the meaning 
of those words of caution which include stress and strain one has to undergo, 
that may lead to a loss of career (academic!). I may repeat the statement of 
Labeyrie [34], 'Recommendation to students: beware of large scale innovative 
projects, they can drive you crazy']. 

In the following, I shall talk here in brief, the basic principle of interferometry 
and its applications, several experiments conducted by the pioneers, as well 
as by us, the problems created by the atmosphere and the remedy, speckle 
interferomctric technique and its use at the telescope in the optical domain, the 
salient features of our newly designed speckle interferometer and the detectors 
are being used for recording speckles etc. The application of non-redundant 



masking technique and its use at the telescopes are also discussed. This part of 
the talk is the concise version of the lecture series which I delivered to the senior 
Ph. D. students at the fall of 1995 at IIA, Bangalore. In addition, I shall discuss 
the various techniques applied to image restoration and their shortcomings, as 
well as the astrophysical problems which can be observed using 2.34 meter VBT, 
at Kavalur using the aforc-mcntioncd interferometer. 

2. Interference of Two Light Waves 

When two light beams from a single source are superposed, the intensity at the 
point of superposition varies from point to point between maxima which exceed 
the sum of the intensities in the beams and minima, which may be zero, known 
as interference [35]. In the following, we shall discuss the degree of correlation 
that exists between the fluctuations in two light waves. 

2.1 The Case of Two Monochromatic Waves: 

The intensity I of light has been defined as the time average of the amount of 
energy which crosses in unit time, a unit area perpendicular to the direction of 
the energy flow. The electric vector E is represented by 



E(r,t)=i?{A(r)e-*'^*} 
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here, r(j. ,^ 2) is a position vector of a point and A is a complex vector aj(r)e'8j(''). 

For a homogeneous plane wave, the amplitudes a'^s are constant, while in phase 
function gj are of the form 

g,(r)=K.r-,5j (2) 

where, K is the propagation vector and (5's are the phase constant, which spec- 
ifies the state of polarization. 
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= - (A^e-^''^* + A*2e2»'^* + 2AA*) (3) 

whence, taking the time average over an interval large compared with the period 
T ^ 2tt/uj 

The intensity I, 



/ =< E^ > 

= ^ (M. p + M. p + M. n 

^-{al + al + al) (4) 

where, < > is the ensemble average and * stands for the complex conjugate. 

Let us suppose now that two monochromatic waves Ei and E2 are superposed 
at a point. The total electric field at that point is 

E = El + E2 (5) 

so that, 

E2 = E^+Ei + 2Ei.E2 (6) 

The intensity at the same point is 

I^Ii+l2 + J12 (7) 

where 

h =< E^ > , I2 =< Bj > 

are the intensities of two waves and 

J12 = 2 < E1.E2 > 
is the interference term. 
Let A and B be the complex amplitudes of the two waves, where 

A, - aie'91 ; B, = he''^' (8) 

The (real) phases Qj and hj of the two waves will be different in general, since the 
two waves will have to travel to the intersecting point by different paths. If the 
same phase difference d is introduced between the corresponding components, 
we have 

gi - hi = g2 - h2 = gs - hs = 6 = —Aip (9) 

Aip is the optical path difference (OPD) between two waves from the common 
source to the intersecting point and Ao is the wavelength in vacuum. 



In terms of A and B 

E1.E2 = - (Ae-*'^* + A*e'"*) . (Be"'"* + B 
1 



- (A.Be-2»'^* + A*.B*e2^'^* + A.B* + A*.B) (10) 



Therefore, 



J12 — 2 < E1.E2 > 

= -(A.B*+A*.B) 

^ {aibi + a2b2 + a^bs) cos S (11) 

This expression shows the dependence of the interference term on the ampUtude 
components and on the phase difference of the two waves. 

If the distribution of intensity resuhing from the superposition of two waves 
propagates in z— direction and hnearly polarized with their E vectors in the 
X— direction, then, 

02 = 03 = ^2 = &3 = 
h = \al ; h = \b\ (12) 



and 



J12 = fli^i cos5 
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The total intensity 



I = Ii+h + 2y/hh cosS (14) 



Irnax = h + h + ^^/hh when \ S \ = 0,2Tr, Att (15) 



Imin = h+ h- 'i-yhh. when \6\ = tt, Stt, Stt (16) 

when /i = I2 

/ = 4/1 cos2- (17) 



The intensity varies between 41 and 0. 

2.2 The Case of Quasi-monochromatic Waves: 

In 1868, Fizeau found the relationship between the aspect of interference fringes 
and the angular size of the light source. The complex degree of co-herence 712 (t) 
of the observed source is defined as follows: 

^"^"^ ^ [rri(0).r22(0)]V2 (18) 

in which 

ri2(r)=<Vi(t).V'2*(i + T)> 

is the inter-correlation function of the field '\\}\ and i/'2, measured at two points 
r\ and r2. 

The ensemble average can be replaced by a time average due to the assumed 
ergodicity of the fields. rii(O) is the average intensity at point r^. 

For quasi-monochromatic sources, the Van Cittert-Zernicke theorem states that 
in quasi-monochromatic light the modulus of the complex degree of coherence of 
the source is equal to the modulus of the normalized spatial Fourier transform 
of the source's brightness [35]. 

, , , I FT of brightness distribution I 
1712(0)1=^ 



total intensity 

0(/l2) I 



0(0) 



(19) 



where, | | stands for the modulus,^ is the Fourier transform and 0{f) is the FT 
of the object at spatial frequency f. 

The modulus of the degree of coherence, sampled at several separations (ri, r2) 
is the visibility function of the source which yields discrete values of the modulus 
of the energy spectrum of the source. 

If both the fields are sent on a quadratic detector, it yields the desired cross-term 
(time average due to time response). The measured intensity at the detector 
would be 

/(ri,r2,r)=<|^i(t)+?A2(t + A) f> 
^h+l2 + 2(/i./2)i/2.i?e{^,2(r)} (20) 



with 

/i=<|^i(t)|2> 

and 

h =<l Mt) ?> 



In order to keep the time correlation close to unity, the delay r must be limited 
to a small fraction of the temporal width Tc- 

AiyTc = l (21) 

where Av is the spectral width and Tc is the temporal width. 

A factor less than unity affects the degree of coherence. The corresponding limit 
for the OPD between two fields is the coherence length, defined by 

Ic = c-Tc = (Ao)VAA (22) 

If T << Tc, we can write, 

7i2(T)=7i2(0)e-2--o- (23) 

where, [vo is mean frequency] 
Let $12 be the argument of 712(7") we have, 

Ir,.-,,r = /i +/2 + 2(/i,/2)i/2i?e{| 712(0) | e**-e-2--o-} (24) 

The measured intensity at a distance x from the origin (point at zero OPD) on 
a screen at distance x from the aperture is 

"27rd(x) 
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(25) 



where d{x) — B.x/{z,X) is the OPD corresponding to x. B is the distance 
between the two apertures. 

The modulus of the spatial coherence of collected fields at the aperture appears 
through the contrast of the fringes which can be measured by the following 
equation. 

T~ 

= 1 712(0) I ^^1^ (26) 

^1 + -'2 

where C is the visibility. 



In order to get | 712(0) |, the measurements of /i and I2 should be made sepa- 
rately. 

3. A Few Experiments to Measure Stellar Diameter 

In order to produce Young's fringes at the focal plane of the telescope, Fizeau 
[1] had suggested to install a screen with two holes on top of the telescope. 
According to him, these fringes remain visible in presence of seeing, therefore, 
allow measurements of stellar diameters with diffraction limited resolution. Ste- 
fan attempted with 1 meter telescope at Obscrvatoire de Marseille and fringes 
appeared within the common Airy disk of the sub- apertures. But he could not 
notice any significant drop of fringe visibility. Since the maximum achievable 
resolution is limited by the diameter of the telescope, he concluded none of the 
observed stars approached 0.1 arcsec. in angular size. 

About half a century later, Michelson, who had spent some years with Fizeau, 
could measure the diameter of the satellites of Jupiter with Fizeau interferometer 
on top of the Yerkes refractor. Similar interferometer (Fizeau mask) was placed 
on top of 100 inch telescope at Mt. Wilson [36] and the angular separation of 
spectroscopic binary star Capella was measured. We too conducted the same 
experiment at VBO, Kavalur using 1 meter Carl-Zeiss telescope and successfully 
recorded fringes of several bright stars with a 16 mm movie camera giving an 
exposure of 16 msec, per frame [21]. 

To overcome the restrictions of the baseline Michelson [37] constructed his stellar 
interferometer by installing a 7 meter steel beam on top of the telescope. It 
was equipped with 4 flat mirrors to fold the beams in periscopic fashion. The 
supergiant star a — Orionis were resolved with this interferometer [38]. 

In this design, the maximum resolution is limited by the length of the girder 
bearing the collectors. The spatial modulation frequency in the focal plane is 
independent of the distance between the collectors. This feature allows to keep 
the same detection conditions when varying the baseline B. The telescope serves 
as correlator, thus, provides zero OPD. 

This experiment had faced various difficulties in resolving stars. These are 
mainly due to the (i) effect of atmospheric turbulence, (ii) variations of refractive 
index above small sub-apertures of the interferometer causing the interference 
pattern to move as a whole, (iii) 7 meter separation of outer mirrors is insufficient 
to measure the diameter of more stars and (iv) mechanical instability prevents 
controlling large interferometer. 

4. Effect of Atmospheric Turbulence 

When an idealized astrophysical source of monochromatic radiation enters in 



the absence of atmosphere, is known as plane wave having uniform magnitude 
and phase across the telescope aperture. The point spread function (PSF) of 
the telescope is the modulus square of the Fourier transform of the aperture 
function. The resolution at the image plane of the telescope is determined by 
the width of the PSF. 

When a flat wavefront passes down through atmosphere, it suffers a phase fluc- 
tuations and reaches the entrance pupil of the telescope with patches of random 
excursions in phase [39]. Due to the motion and temperature fluctuations in 
the air above the telescope aperture, inhomogencities in the refractive index 
develop. These inhomogencities have the effect of breaking the aperture into 
cells with different values of refractive index that are moved by the wind across 
the telescope aperture. Kolmogrov law represents the distribution of turbule 
sizes, from millimeters to meters, with lifetimes varying from msecs to seconds. 
Changes in the refractive index in different portions of the aperture result to the 
phase changes in the value of the aperture function. The time evolution of the 
aperture function implies that the PSF is time dependent. If the atmosphere 
is frozen at a particular instant; each patch of the wavefront with diameter Tq 
— Fried parameter — would act independently of the rest of the wavefront re- 
sulting in many bright spots. These spots are known as speckles and spread 
over the area defined by the long exposure image. Computer simulated analysis 
demonstrates the destructions of the finer details of an image of a star by the 
atmospheric turbulence [23] . The size of To is found to be varied between 8 to 
12 cm at Hq, wavelength during the night at the 2.34 meter VBT [40]. 

The long-exposure PSF is defined by the ensemble average, < S'(x) >, indepen- 
dent of any direction. The average illumination, /(x) of a resolved object, 0(x) 
obeys convolution relationship, 

< /(x) >= ©(x)* < 5(x) > (27) 

where, * stands for convolution, x = (x,y) is a 2-dimcnsional space vector. Using 
2-dimcnsional Fourier transform, this equation can be read as, 

< /(u) >= 6(u)- < 5(u) > (28) 

where, 0(u) is the object spectrum, < S{u) > is the transfer function for long- 
exposure images and is the product of the transfer function of the atmosphere 
B{u), as well as the transfer function of the telescope, T{u.). u is the spatial 
frequency vector with magnitude u. The transfer function for long-exposure 
image can be expressed as, 

< S{u) >= B{u) ■ r(u) (29) 

The benefit of the short-exposure images over long-exposure can be visualized by 
the following explanation. Let us consider two seeing cells separated by a vector 
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in the telescope pupil, Au, where A is the mean wavelength and u is an angular 
spatial frequency vector. If a point source is imaged through the telescope 
by using pupil function consisting of two apertures, corresponding to the two 
seeing cells, then a fringe pattern is produced with narrow spatial frequency 
bandwidth. If the major component /(u) at the frequency u is produced by 
contributions from all pairs of points with separations Au, with one point in 
each aperture and is averaged over many frames, then the result for frequencies 
greater than Tq/A tends to zero. The Fourier component performs a random 
walk in the complex plane and average to zero, < /(u) > = 0, when u > Tq/A. 

For a large telescope, the aperture, P, can be sub-divided into a set of sub- 
apertures. Pi. According to the diffraction theory [35], the image at the focal 
plane of the telescope is obtained by adding all such fringe patterns produced 
by all possible pairs of sub-apertures. With increasing distance of the baseline 
between two sub-apertures, the fringes move with an increasingly larger ampli- 
tude. On a long-exposure images, no such shift is observed, which implies the 
loss of high frequency components of the image. While, in the short-exposure 
images (< 20msec), the interference fringes are preserved. 

5. Speckle 

The term 'Speckle' refers to a grainy structure observed when an uneven surface 
of an object is illuminated by a fairly coherent source. A good example of 
speckle phenomena may be observed at the river port when many boats are 
approaching towards the former at a particular time or in the swimming pool 
when many swimmers are present. Each boat or swimmer emits wave trains and 
interference between these random trains causes a speckled wave field on the 
water surface. Depending on the randomness of the source, spatial or temporal, 
speckles tend to appear. Spatial speckles may be observed when all parts of 
the source vibrate at same constant frequency but with different amplitude 
and phase, while temporal speckles are produced if all parts of it have uniform 
amplitude and phase. With a non-monochromatic vibration spectrum, in the 
case of random sources of light, spatio-temporal speckles are produced. 

The ground illumination produced by any star has fluctuating speckles, known 
as star speckles. It is too fast and faint, therefore, cannot be seen directly. 
Atmospheric speckles can be observed easily in a star image at the focus of 
a large telescope using a strong eyepiece. The star image looks like a pan of 
boiling water. If a short exposure image is taken, speckles can be recorded. The 
speckle size is of the same order of magnitude as the Airy disc of the telescope 
in the absence of turbulence. The number of correlation cells is determined by 
the equation N — D/tq- As the seeing improves, the number decreases. 
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5.1 Imaging in the Presence of Atmosphere : 

Let us consider an imaging system consists of a simple lens based telescope in 
which the point spread function (PSF) is invariant to spatial shifts. An object 
(point source) at a point x' anywhere in the field of view will, therefore, produce 
a pattern S'(x — x') across the image. If the object can emit incoherently, the 
image /(x) of a resolved object 0(x) obeys convolution relationships. The 
mathematical description of the convolution of two functions is of the form: 

/(x) = / 0(x)'S'(x - x')dx' (30) 

Convolution equations can be reduced to agreeable form using the Fourier con- 
volution theorem. The Fourier transform of a convolution of two functions is the 
product of the Fourier transform of the two functions. Therefore, in the Fourier 
plane the effect becomes a multiplication, point by point, of the transform of 
the object 0(u) with the transfer function S'(u). 

5.2 Laboratory Simulation : 

Atmospheric seeing can be simulated at the laboratory by introducing distur- 
bances in the form of a glass plate with silicone oil [2]. We had introduced 
various static dielectric cells (SDC) of various sizes etched in glass plate with 
hydrofluoric acid. Several glass plates with both regular and random distribu- 
tion of SDCs of known sizes were made and used in the experiment [31]. The 
phase-differences due to etching lie between 0.2A and 0.7A. In order to obtain 
the light beam from a point source, similar to the star in the sky, we had devel- 
oped an artificial star image by placing a pair of condensing lenses along with 
micron-sized pin-hole in front of the source [22] . The beam was collimated with 
a good quality Nikkon lens; the wave fronts from this artificial star enter a sim- 
ulated telescope whose focal ratio is 1:3.25 (similar to the prime focus of VBT). 
The image was magnified to discern the individual speckles with a high power 
microscope objective. The speckles were recorded through a lOnm interference 
filter centered on 5577A. Figure 1 depicts (a) the laboratory set up to simulate 
speckles from an artificial star, (b) speckles obtained in the laboratory through 
the aforementioned narrow band filter. The image was digitized with the PDS 
lOlOM micro-densitometer and processed using the COMTAL image processing 
system of the VAX 11/780 at the VBO, Kavalur. The clipping technique was 
used to enhance the contrast in grey levels. The clipped image is superposed 
on the histogram-equalized original image. The laboratory set up was sensi- 
tive enough to detect aberrations produced by the objective lenses, as well as 
micro-fluctuations in the speckle pattern caused by vibrations. By introducing 
an aperture mask in front of this telescope, we could obtain the fringes [22, 
29, 30]. The similarity of the observed image shape at the laboratory with the 
computer simulations was found. 
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Figure 1: (a). Laboratory set up to simulate speckles from an artificial star, 



The laboratory simulation is necessary for the accurate evaluation of the 
performance of speckle imaging system; comparison of the experimental results 
should be made with the computer simulations. The importance of the system- 
atic use of simulated image is to validate the image processing algorithms in 
retrieving the diffraction limited information. 

5.3 Speckle Interferometer : 

A speckle interferometer is a high quality diffraction limited camera where mag- 
nified (~ f/100) short exposure images can be recorded. Additional element for 
atmospheric dispersion corrections is necessary to be incorporated. At an in- 
creasing zenith distance speckles get elongated owing to this effect. Either a pair 
of Risley prism must be provided for the corrections or the observation may be 
carried out using a narrow bandwidth filter. In the following, the salient features 
of our newly developed speckle interferometer [27, 28] are described in brief. 

The wave front falls on the focal plane of an optical flat made of low expansion 
glass with a high precision hole of aperture (~350 /i), at an angle of 15° on 
its surface [27]. The image of the object passes on to the microscope objective 
through this aperture, which slows down the image scale of this telescope to 
f/130. A narrow band filter would be placed before the detector, to avoid the 
chromatic blurring. The surrounding star field of diameter 10 mm, gets reflected 
from the optical flat on to a plane mirror and is re- imaged on an intensified CCD, 
henceforth ICCD [24]. We have recorded a large numbers of specklegrams of 
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several close binary systems and of other point source using uncooled ICCD as 
sensor. The image at the Cassegrain focus of 2.34 meter VBT is sampled to 0.015 
arcsec per pixel. Figure 2 depicts (a) the optical layout of the interferometer, 
(b) speckles of the close binary star HR4689 obtained with the same set up at 
Cassegrain focus of the said telescope on 28*''-l^* March, 1997. 



This interferometer is built with extreme care so to avoid flexure problems 
which might affect high precision measurements of close binary star systems 
etc., in an unfavourable manner. The design analysis has been carried out with 
the modern finite element method [41] and computer aided machines were used 
in manufacturing to get dimensional and geometrical accuracies. The method 
requires the structure to be subdivided into a number of basic elements like 
beams, quadrilateral and solid prismatic elements etc. A complete structure 
is built up by the connection of such finite elements to one another at a def- 
inite number of boundary points called nodes and then inputting appropriate 
boundary constraints, material properties and external forces. The relationship 
between the required deformations of the structure and the known external 
forces is [i4r]{d} — {F}, where, \K\ is the stiffness matrix of the structure, {d} 
is the unknown displacement vector and {F} is the known force vector. All the 
geometry and topology of the structure, material properties and boundary con- 
ditions go into computation of \K\. The single reason for universal application 
of finite clement method is the ease with which the matrix, \K\ is formulated 
for any given structure. 
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Figure 2: (a). Optical layout of the speckle interferometer, (b). Speckles of the 
close binary star HR4689 obtained with the new speckle interferometer at VBT. 
The numbers on the axes denote pixel numbers with each pixel being equal to 
0.015 arcsec. 
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5.4 Aperture Synthesis : 

Two promising methods, viz., (a) speckle masking (see chapter 7.3) technique 
[42], (b) non- redundant aperture masking technique [6, 7] are based on the 
principle of phase-closure method. To recover the Fourier phases of the source 
brightness distribution from the observations, it is necessary to detect fringes on 
a large baselines, therefore, enables one to reconstruct images. The concept of 
using three antennae arranged in a triangle was first introduced in radio astron- 
omy [43] in late fifties. Closure-phases are insensitive to the atmospherically 
induced random phase errors, as well as to the permanent phase errors intro- 
duced by the telescope aberrations in optics. Since any linear phase term in the 
object cancels out [44], this method is insensitive to the position of the object 
but sensitive to any object phase non-linearity. 

The measurements of the closure-phases was first obtained at high light level 
with three- hole aperture mask placed in the pupil plane of the telescope [5]. 
Interference patterns of the star were recorded using CCD as sensor. We had 
conducted similar experiment around the same time by placing an aperture 
mask of 3-holcs, 10cm in diameter, arranged in a triangle, over 1 meter telescope 
at VBO, Kavalur and tried to record the interference pattern with our earlier 
version of the interferometer [21]. In this experiment, we used a 16mm movie 
camera as detector. The advantage of placing aperture mask over the telescope, 
in lieu of pupil mask is to avoid additional optics. But a curious modulation of 
intensity in the fringe pattern was noticed, therefore, unable to proceed further. 
The modulation could result from a time-independent aberration in the optical 
system of the telescope [21]. However, we performed the similar experiment in 
the laboratory, using the aperture mask of 3-holes, placed in front of a simulated 
telescope and were able to record the interference patterns [30]. 

The aperture synthesis imaging technique with telescope involves observing an 
object through a masked aperture of several holes and recording the interference 
patterns in a series of short-exposure. The patterns contain information about 
structure of the object at the spatial frequencies from which an image of the 
same can be reconstructed by measuring the visibility amplitudes and closure 
phases. This method produces images of high dynamic range, but restricts to 
bright objects. 

Several groups have obtained the fringe patterns using both non-redundant and 
partially redundant aperture mask of N-holes at large or moderate telescope 
[8-15]. In the laboratory, the shapes of fringe pattern of N-hole apertures were 
also studied by us by introducing the various aperture masks arranged in both 
non-redundant and partially redundantly [22, 29]. A few instruments developed 
for the pupil plane aperture mask are described below. 

(i) The instrument developed for the Hale 5 meter telescope [11] used f/2.8 {4> 
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85inni) Nikkon camera lens to coUimate f/3.3 primary beam of the telescope. 
The lens formed an image of the primary mirror at a distance of about 85 mm 
where a mask was placed on a stepper-motor-driven rotary stage controlled by 
a PC. Another identical lens forms a second focus (scale is 12 arcsec/mm). This 
image was expanded by microscope objective (X80), enabhng to sample 0.15 
arcsec/mm on the detector. A narrow band interference filter (6300A, FWHM 
30A) was placed between the microscope objective and the detector. The re- 
sistive anode position sensing photon counting detector [45] was used to record 
the interference patterns. They were successful in producing optical aperture 
synthesis maps of two binary stars [3 — CoronaBorealis and a — Hercules. 

(ii) University of Sydney [12] had developed a masked aperture-plane inter- 
ference telescope (MAPPIT) for the 3.9 meter Anglo-Australian telescope to 
investigate interferometry with non-redundant masks. A field lens re-images 
the telescope pupil down to diameter of 25 mm and the aperture mask is placed 
where the pupil image is formed. Dove prism is used to rotate the field, allowing 
coverage of all position angles on the sky. Dispersed fringes are produced using 
a combination of image and pupil plane imaging. The camera lens and micro- 
scope objective produce an image of a star in one direction. In the orthogonal 
direction, the detector receives the dispersed pupil image. The mask holes play 
the role of the spectrograph slit. A cylindrical lens is used as the spectrograph's 
camera lens. Image photon counting system is used in this experiment to record 
the dispersed fringe pattern. This instrument is used at coude focus of 3.9 meter 
Anglo-Australian telescope. They were able to resolve several close binaries, as 
well as to measure angular diameters of cool stars [13, 14]. 

Recently, Bedding [15] of the said University has developed another version of 
this technique, called multiplexed one-dimensional speckle [MODS], by replacing 
holes in the aperture mask with a slits. Using a cylindrical lens that creates a 
continuous series of one-dimensional intcrfcrograms, intcrfcrograms from many 
arrays can be recorded side- by-side on a 2-dimcnsional detector. A narrow 
band filter is used in place of dispersing prism. The mask containing slit is 
placed in the coUimating beam. The optics in the interference direction form 
a image plane interferometer and in the orthogonal direction, the cylindrical 
lens produces a pupil image. Measurements at different position angles can 
be made by rotating this lens. Two-dimensional image reconstruction can be 
made from the series of one-dimensional intcrfcrograms [46]. Since the same 
number of photons will be collected in this technique as the conventional speckle 
interferometer, the limiting magnitude would not change [15]. 



5.5 Detector : 

The speckle interferometer requires snap shots of very high time resolution of 



19 



the order of (a) frame integration of 50 Hz [47], (b) photon recording rates of 
a few MHz [48]. We have been using an uncooled ICCD as sensor to record 
speckle-grams of the objects. It gives video signal as output. The images were 
recorded at an exposure 20 msec using a frame grabber card DT-2861 manu- 
factured by Data Translation'^^. Each frame consists of odd and even field. 
But the coherence time of the atmosphere is a highly variable parameter [49] , 
it is desirable to observe speckles with a photon counting system. In the follow- 
ing, we shall describe the salient features of precision analogue photon address 
(PAPA) detector. 

The PAPA is a 2-d photon counting camera [48]. Individual photons hitting 
the photocathode of a high gain image intensifier produce a spot of light from a 
phosphor screen on the output side of the intensifier. An image of the phosphor 
screen is sent through optics to 19 photo- multiplier tubes (PMT), 18 of which 
have their active area covered by one of 9 different grey scale masks. The 
19th tube acts as an event strobe, registering a digital pulse if the spot on the 
phosphor is detected by the instrument. 9 tubes are used to obtain positional 
information for an event in one direction, while the other 9 are used for positional 
information in the orthogonal direction. If the phosphor hit in a region where it 
is not covered by the mask, an event is registered by the phototube. By taking 
the information from each phototube for each event, these are recorded as a list 
of photon addresses and arrival times in a binary format. 

6. Data Processing 

Stellar speckle interferometry consists of taking many short integration images 
of an object. The frame image is the convolution of the instantaneous PSF 
of telescope— atmosphere 5(x) with the actual object irradiance distribution 
0(x). The intensity distribution /(x) of the speckle interferograms, in the case 
of quasi-monochromatic incoherent source can be described by the following 
space- invariant imaging equations. 

/(x) = 0(x) * S(x) (31) 

In the Fourier domain, convolution becomes an ordinary product so that, 

7(u) ^ 0(u) • S(u) (32) 

The ensemble averaged power spectrum is given by, 

<|/(u)|^>=|6(u)|^.<|S(u)|^> (33) 

Since | 5'(u) |^ is a random function in which the detail is continuously chang- 
ing, the ensemble average of this term becomes smoother. The smooth function 
can be performed on a point source yields <| 5(u) |^>. The object autocor- 
relation can be obtained by inverse Fourier transform. Figure 3 depicts the 
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Figure 3: Auto-correlation of a binary star HR5747. The speckles were recorded 
with earlier version of the speckle interferometer [21, 24] on 16/17th March 1990. 
The numbers on the axes denote pixel numbers with each pixel being equal to 
0.026 arcsec. 



auto-correlation of a binary star HR5747. The specklegrams of this star were 
obtained at the Cassegrain focus of the 2.34m VBT using speckle interferometer 
[21] on 16/17th March 1990. The uncooled ICCD was used as sensor to record 
the specklegrams; each pixel of this ICCD was sampled to 0.026 arcsec. 



7. Image Processing 

Several methods have been developed to produce a map of the phase portion 
of the diffraction limited object Fourier transform. In the following we shall 
discuss some of the techniques. 

7.1 Speckle Holography : 

I. In speckle holographic technique [50] complete image reconstruction becomes 
possible when a reference print source is available in the field of view (within 
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isoplanatic patch < 7 arcsec). 

Let us represent the point source by a Dirac impulse AS{x.) at the origin and 
let Oi(x) be a nearby object to be reconstructed. The irradiance distribution 
in the field of view is 

0(x)=A<5(x) + Oi(x) (34) 

A regular speckle intcrferometric measurement will give the squared modulus of 
its Fourier transform 0(u) 

|6(u)p=.|A + 6i(u)r 

= A^ + Adi (u) + a6* (u) + 61 (u)6* (u) (35) 

The inverse Fourier transform gives the autocorrelation Co{^) of the field of 
view 

Co(x) = A^Si^) + ^Oi(+x) + AOi(-x) + Coi(x) (36) 

where Coi(x) is the autocorrelation of the object. 

The first and the last term in equation (36) arc centered at the origin. If the 
object is far enough from the reference source, 0(x) its mirror image 0(— x) is 
therefore, recovered apart from a 180° rotation ambiguity. 

II. The alternate idea is to calculate the average cross spectrum between the 
objects and the reference. Let Oi(x) and 02(x) be respectively the object and 
reference brightness distributions and /i (x) and I2 (x) their associated instanta- 
neous image intensity distributions. The relation between the objects and the 
images in the Fourier space becomes 

/i(u) = 6i(u).S(u) (37) 

/2(u) = 62(u).S(u) (38) 

where, 5(u) is the impulse response of the telescope and the atmosphere. 
The average cross-spectrum between the object and the reference 

< /i(u)/*(u) >- di{u)d;{u)- <\ 5(u) p>M (39) 

The equation (39) shows that as the speckle holography transfer function <| 
S'(u) p> is real, the method is insensitive to aberrations and the phase of the 
cross spectrum expected value coincides with the phase difference between the 
object and the reference. 
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7.2 Knox-Thomson Technique (KT) : 

The method [51] involves in interfering /(u) with itself after translating by a 
small shift vector Au. The KT correlation may be defined in Fourier space as 
products of /(u). 

/(u, Au) =< /(u)r (u + Au) > (40) 

This gives us the product, 

</(u)/*(u + Au) >=6(u)6*(u + Au)-<§(u)S^*(u + Au) > (41) 

Introducing the phases explicitly in the form O =| O | e^'^° etc. and using 
A(/) ~ (piu) — (/'(ui + Au), we have 

/(u) II /(u+Au) I e^''" ^\ 6(u) II 0(u+Au) | e^'^"- | S{u) \\ S{u+Au) \ e^'"^^ 

(42) 
If this equation is averaged over a large number of frames, the feature (Ai^s) = 
0. When Au is small, | 0(u + Au) | « | 0(u) | etc. and so 

< /(u)r (u + Au) >=| 6(u) p e^'^'^o • (a smooth function) (43) 
from which, together with equation (32), Ai/iq can be determined. 

7.3 Speckle Masking or Triple Correlation Technique (TC) : 

In case of non-availability of reference point source within iso-planatic patch, 
the instantaneous PSF can be estimated from the speckle pattern itself. Weigelt 
[52] suggested to multiply the object speckle pattern /(x) by an appropriately 
shifted version of this /(x + Xi). For a binary star, the shift is equal to the 
angular separation between the stars, masking one of the two component of 
each double speckle. The result is correlated with /(x). The Fourier transform 
of the triple correlation is called bispectrum and its ensemble average [42] is 
given by 

/(U1,U2) =< /(Ui)r(ui +U2)/(U2) > (44) 

where, /(ui), /(U2), /(ui + U2), denote the Fourier transforms of /(x). 

In the second order moment or in the energy spectrum, phase information is 
lost, but in the third order moment or in the bispectrum it is preserved. If we 
put equation (32) into equation (44), it emerges as, 

/(Ui, U2) - 0(ui)6*(ui + U2)6(U2) < 5(ui)5*(ui + U2)5(U2) > (45) 
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The relationship shows that the image bispectrum is equal to the object bis- 
pectrum times a bispectral transfer function. The object bispectrum is given 

by, 



/o(ui,U2)-0(ui)0*(ui+U2)0(u2) 

^ < J(Ui)7*(Ui+U2)/(u2)> 
<Siui)S*{ui+U2)S{u2)> 



(46) 



The modulus | 0(u) | of the object Fourier transform 0(u) can be derived from 
the object bispectrum /o(ui,U2) [53]. 

The phase 0(u) of the object Fourier transform can also be derived from the 
object bispectrum. Let (j)b be the phase of the object bispectrum and we get, 

6(u) =1 6(u) I e^'^(") (47) 

and 

To{u,,U2) =1 /o(ui, U2) I e-''*''("i'"^) (48) 

Equations (47) and (48) may be inserted into equation (42), therefore, yields 
the relations, 



7o(ui,U2) =1 6(ui) I e^-^f^i) I 6(U2) I e^-^f"^) | 6(ui + U2) | e-^'^("i+"^) ^ 

(49) 

gj06(ui,U2) _ gi0(ui)gj0(u2)g-j0(ui+U2) _^ /gQ\ 

0f,(ui, U2) = (f>{ui) + (f){u2) - (f>{ui + U2) -^ (51) 

(j>{ui + U2) = (/)(u) 
= <?!)(ui)+<?i(u2)-</'6(ui,U2) (52) 

Equation (52) is a recursive equation for calculating the phase of the object 
Fourier transform at coordinate u = Ui + U2 [53]. If the object spectrum at Ui 
and U2 are known, the phase spectrum at (ui + U2) can be computed. But the 
bispectrum phases are of mod 27r, therefore, the reconstruction in equation (48) 
may lead to tt phase mismatches between the computed phase-spectrum values 
along the different paths to the same point in frequency space. The unit am- 
plitude phasor recursive re-constructor are insensitive to the phase ambiguities 
and the computing argument of the term, e^'^^"'^^^^' , can be expressed as, 
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gJ0(ui+U2) _ gj[i*(ui) + 0(u2)-06(ui,U2)] ('53') 

The triple correlation algorithm based on this unit amplitude recursive re- 
constructor method was developed at our institute for processing the stellar 
objects [33], as well as for the extended object [54]. 

7.4 Relationships : 

From the sections 7.2 and 7.3, we find the relationship among the two widely 
used algorithms, namely, KT and TC methods [55]. In autocorrelation tech- 
nique, the major Fourier component of the fringe pattern is averaged as a prod- 
uct with its complex conjugate and so the atmospheric phase contribution is 
eliminated and the averaged signal is non-zero (see chapter 4). Unfortunately 
the phase information is not preserved. 

The KT is a small modification of the autocorrelation technique In the KT 
method, approximate phase closure is achieved by two vectors, u and u + Au, 
and assuming that the pupil phase is constant over Au. The major Fourier 
component of the fringe pattern is averaged with a component at a frequency 
displaced by a vector Au. The vector displacement Au should not force the 
vector difference — u — Au outside the spatial frequency bandwidth of the fringe 
pattern which, in turn, preserves the Fourier phase difference information in the 
averaged signal. The atmospheric phase effectively forms a closed loop. 

In the bispectrum method a third vector, Au is added to form phase closure. 
When AAu > Tq, the third vector Au is essential. The KT method fails with 
this arrangement. In this system, the two apertures are extended to three and 
the Au > To/A. It can be seen that atmospheric phase contribution is not closed. 
KT is limited to frequency differences Au < r^/A. If the bispectrum average 
is performed, the phase is closed and Fourier phase difference information is 
preserved. This method can obtain phase information for phase difference Au > 
ro/A. 

7.5 Blind Iterative Deconvolution Technique (BID) : 

In this technique, the iterative loop is repeated enforcing image-domain and 
Fourier-domain constraints until two images are found that produce the input 
image when convolved together [56, 57]. The image-domain constraints of non- 
negativity is generally used in iterative algorithms associated with optical pro- 
cessing to find effective supports of the object and or PSF from a speckle-gram. 
The implementation of the BID, used by us is described below. 

The algorithm has the degraded image /(x) as the operand. An initial estimate 
of the PSF, S'(x), has to be provided. The degraded image is deconvolved 
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from the guess PSF by Wiener filtering, which is an operation of multiplying 
a suitable Wiener filter ( constructed from the Fourier transform S'(u) of the 
PSF) with the Fourier transform J(u) of the degraded image. The technique 
of Wiener filtering damps the high frequencies and minimizes the mean square 
error between each estimate and the true spectrum. This filtered deconvolution 
takes the form 

6(u)=/(u)^ (54) 

5(u) 

The Wiener filter, 0/(u), is given by the following equation: 

a,(u) . . ^(")^*(-^) (55) 

I 5(u) P + I iV(u) P 

This noise term, N{u), can easily be replaced with a constant estimated as the 
rms fluctuation of the high frequency region in the spectrum, where the object 
power is negligible. The Wiener filtering spectrum, 0(u), takes the form: 

0(u) = I{u)^ ^ ^^ ^ (56) 

S{u)S*{u)+N{u)N*{u) 

This result 0(u) is transformed to image space, the negatives in the image are set 
to zero, and the positives outside a prescribed domain (called object support) 
are set to zero. The average of negative intensities within the support are 
subtracted from all pixels. The process is repeated until the negative intensities 
decreases below the noise. 

A new estimate of the PSF is next obtained by Wiener filtering the original 
image /(x) with a filter constructed from the constrained object 0(x). This 
completes one iteration. This entire process is repeated until the derived values 
of 0(x) and S{x) converge to sensible solutions. Before applying this scheme 
of BID to the real data, we have tested the algorithm with computer simulated 
convolved functions of binary star and the PSF caused by the atmosphere and 
the telescope. The reconstruction of the Fourier phase of these are shown in 
figure 4. The starting guess for the PSF used for this calculation was a Gaussian 
with random noise. We were able to obtain the output image, as well as the 
output PSF after 225 iterations. 



The uniqueness and convergence properties of the deconvolution algorithm are 
uncertain for the evaluation of the reconstructed images if one uses the BID 
method directly. We have tested the code [58] and found that it is essential to 
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Figure 4: 2-dimcnsional maps of the simulated image of (a) binary star, (b) 
simulated atmospheric PSF, (c) the convolved functions of these two (bpsf), (d) 
the retrieved image of the binary star, (e) reconstructed atmospheric PSF. 
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estimate the input support radius of the object using the auto-correlation tech- 
nique, which helps in completing the reconstruction satisfactorily after several 
iterations [32]. 

8. Astrophysical Programmes 

The contribution of single aperture interferometry to several important fields in 
astrophysics has increased considerably, viz., separation and orientation of close 
binary stars [32, 59-61], shapes of asteroids [62], mapping of the finer features 
of extended objects [26], sizes of certain types of circumstellar envelopes [63, 64], 
structures of active galactic nuclei [65, 66], resolving the gravitationally lensed 
QSO's [67] etc. We too plan to observe the following interesting objects, if an 
intensified photon counting detector [48] can be made available. 

(i) Studies of binary stars play a fundamental role in measuring stellar masses, 
providing the benchmark for stellar evolution calculations. A long term benefit 
of speckle imaging is a better calibration of the main-sequence mass-luminosity 
relationship. Most measurements have been made at large aperture telescopes 
by groups in France, Russia, and the United States. Programmes of binary 
star interferometry are being carried out at telescopes of moderate and small 
aperture too. But measurements from the southern hemisphere continue to be 
rare. Many rapidly moving southern binaries are being ignored and a number 
of discoveries are yet to be confirmed. 

(ii) Most of the late-type stars are available in the vicinity of sun. All known 
stars, within 5 pc radius from the sun are red dwarfs with m^ > +15. Due to 
the intrinsically faint nature of K- and M- dwarfs, their physical properties are 
not studied extensively. These dwarfs may often be close binaries which can be 
detected by speckle interferometric technique. 

(iii) Another important field of observational astronomy is the studies of the 
physical processes, viz., temperature, density and velocity of gas in the active 
region of the active galactic nuclei ( AGN) . Optical imaging in the light of emis- 
sion lines on sub-arcsec scales can reveal the structure of the narrow-line gas. 
The scale of narrow-line regions is well resolved by the diffraction limit of a 
moderate-sized telescope. 

(iv) The spatial distribution of circumstellar matter surrounding objects which 
eject mass, particularly young compact planetary nebulae can be determined 
[63, 64]. 

(v) Capability of resolving the gravitationally lensed QSO's in the range of 0.2 
arcsec. to 0.6 arcsec. will allow the discovery of many more lensing events [67]. 
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9. Epilogue 

The understanding of the basic random interference phenomenon — speckle — 
is of paramount importance to the observational astronomy. In recent years 
the uses of speckle pattern and a wide variety of applications have been found 
in many other branches of physics and engineering. Though the statistical 
properties of speckle pattern is complicated, detailed analysis of this is useful in 
information processing. Though the stellar speckle interferometry is capable of 
detecting relatively faint objects (~ 16th. magnitude), the angular resolution 
is limited by the diameter of the telescope. Angular resolution can vastly be 
improved by using long base line interferometry using two or more telescopes. 

India has an outstanding group in the field of radio astronomy using long base- 
line interferometry. From the experience we have gained in developing the field 
of optical interferometry, we are confident in building a long baseline interfer- 
ometer in the optical and IR band. Long baseline interferometric observations 
of the objects [34, 68] would offer the possibilities for direct measurement of all 
the basic physical parameters for a large number of stars. 
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